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ABSTRACT 

The problem of one pair of identical nucleons sitting in M 
single particle levels of a potential well and interacting through 
the pairing force is treated introducing, in the Hamiltonian for- 
malism, even Grassmann variables. The eigenvectors are analyt- 
ically expressed solely in terms of these with coefficients fixed by 
the eigenvalues and the single particle energies. For these a spe- 
cific model is needed: in the case of the harmonic oscillator well, 
for any strength of the pairing interaction, an accurate expression 
is derived for both the collective eigenvalue and for those trapped 
in between the single particle levels. Notably the latter are labelled 
through an index upon which they depend parabolically. 
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We have recently obtained, in the framework of the Grassmann algebra, 
the analytic expressions for the eigenvalues and the eigenvectors of n pairs of 
like-nucleons interacting through the pairing Hamiltonian and sitting in one 
single-particle level 

When the single-particle levels, with angular momentum • • • jjv and 
energies e±, e 2 , ■ • • ejy (all the j's being assumed to be different), are M usually 
the problem is dealt with numerically. In this case the Hamiltonian of the 
system reads 



JV j v JV in jv 

H = E e " E ^j v m v \ v m v ~ G ^2 ^j^m^^j^rhf, E \ v m v \ v m v i 

u=l m„=—j v /i,i/=l m M =l/2 m„=l/2 

(1) 

where Xj m and X* m are the odd (anticommuting, nilpotent) Grassmann vari- 
ables (associated to the nucleons annihilation and creation, respectively) and 



Xjm = (— 1) J m Xj- m ■ (2) 

Introducing even (commuting, nilpotent) Grassmann variables to describe 
a pair of fermions with vanishing third component of the total angular mo- 
mentum (M=0), namely || 

'••pjrn = ( 1 )^ Xj— m \j m , (3) 

(HD becomes 



H =H e v E X )vmvKrnv - G E E f^m.fjvruv ■ (4) 

u=l m v =-j v n,c=l m )1 =l/2m„=l/2 

In this letter we provide the eigenvectors of one pair of fermions interact- 
ing through the Hamiltonian (|4]) in the presence of M single particle levels 
in terms of ([J) with coefficients expressed through the eigenvalues and the 
single particle energies. When the latter are those of an harmonic oscillator, 
an accurate, almost analytic expression is given for all the eigenvalues as 
well. 

We start by counting the total number N s of states available to the sys- 
tem: it is given by 
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where 

^ = £^ = £(^ + 1/2) • (6) 

u=l u=\ 

In the above the states with the two fermions sitting on two different levels 
are those with the two fermions placed on the same level are . 

Since we are interested in the physics where the pairing force is active, we 
consider of the states only the ones having the two fermions in time 
reversal orbits: these are Y^f=i 

Notwithstanding the presence of both the A's and the <^'s in (f|), we search 
for eigenstates in the form 

JV 3v 

^ = £ £ fij v rn v tp) vmv ■ (?) 
v=\ m v =\j1 

Using the Grassmann algebra rules, the coefficients (3 are then found to obey 
the system of fl equations 



(£-2e u )f3 Jvmv + J2 £ 0^ = 0, (8) 

M=l m„=l/2 

where 1 < v < Af, 1/2 < m„ < j u . 

We cast the above system in the Q x Q matrix form 



/ fill 


Byi ■ 


■ B w \ 


/ Pjil/2 


\ 


B21 


B22 ■ 


• B 2 x 


Pji3/2 






BM2 ■ 


■ -Saw / 


V 03N3H 


) 



where the elements of the blocks B^ u , of dimensions Q u x Q^, are [B^ v ] a p = 
(S - 2e v ) 5^ u S Q/ 3 + 1. 
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In the above all the energies are expressed in units of G, i.e. 



e„ = e u /G and £ = E/G . (10) 
The matrix M is easily diagonalized, the associated determinant being 



dct M 



P =i 



2e, 



Af Af Af 

n(£-2<g+£ n „ n (£- 2 ^ 



'ID 



and the corresponding secular equation admits 2J\f distinct solutions. Clearly 
M of these are 



£ p = 2e p , with degeneracy d p = Q p — 1 , 1 < p < J\f , 
the remaining Af (non-degenerate) fulfilling instead the equation || 



Af 



l + f(£)=0 with /(£) = £ 



1 f - 2e„ 



;i2) 



(13) 



The eigenvalues (|12D clearly correspond to states insensitive to the pairing 
interaction, being associated to a pair with non vanishing angular momen- 
tumF]. The solutions of (|T3"D correspond instead to states (describing a pair 
coupled to an angular momentum J = 0) affected by the pairing force. 

Now, beyond the trivial solutions (0), also the eigenvalues fulfilling (|13|) 
and the associated eigenstates can be analytically obtained (at least, as we 
shall see, for an harmonic oscillator potential well) in a basis, of dimensions 
lower than (^), which generalizes the one we introduced in ref.JT[]. 

To show how this occurs we organize, as in ref.0, the states associated 
with each of the Af single particle levels into two sets, embodying Q v — 1 
and 1 levels, respectively (of course it must be Q„ ^ 1). Correspondingly we 
define the following new basis of 2jV normalized states with 1 < v < jV 
and a = 0, 1: 



$(o) 

v 

$(1) 



_J_y^>-l , n . 



V 



■N 



(14) 



* Obviously if f2„ = 1, then d v = 0: hence the free solution ( |12| ) is absent. Indeed a 
pair of fermions on the level j u = 1/2 can only couple to angular momentum J = 0, hence 
feeling the pairing interaction. 
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In this basis the rectangular blocks B pv become 



B 



£>(0) 



a 



(15) 



where the row vector a M , of dimension (f2^ — 1), is filled with ones and Bj$ 
and Bj$ are two matrices of dimensions (0„ — 1) x (f^ — 1) and lxl, 
respectively. 

Then the 27V x 2JV matrix Ai, representing the operator 7i = £ — H/G 
in the basis (0), is obtained by replacing each B^ with a 2 x 2 matrix with 
elements given by the sum of the elements of each block in ( |T5| ) (divided by 
the corresponding normalization factors). One thus gets 



( £ -2e l + d l 

£ - 2ei + 1 



M 



\fd~l 



whose determinant reads 



1 



\/d x d M 



£ - 2e^ + djv 



\ 



ydM 
£ - 2e N + 1 / 



(16) 



detM 



det 



V 



1 + £ - 2ei 



AT 

11(^-2^ 



1 , £-2eAf 



i+£-2^ ; 



AT AT AT 

n(£-2e M )+£^ n (^-2^ 

1-1=1 V=\ ^1/=1 



AT 



(17) 



From C pTrl) both the Af unperturbed solutions £ p = 2e p and the TV ones 
affected by the pairing interaction follow. 

The latter are also obtained in the normalized basis (introduced long ago 
by Richardson H with a different technique) 
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$(°) + 



v 



(18) 



which linearly combines the building blocks of (|14l) . Note that the index of 
nilpotency of the commuting variables Q u is Q u : hence the (|1^) are "more 
bosonic" than the (§). Indeed they are often referred to as s-quasibosons 
( J = 0) and lead to the following TV dimensional representation of 7i f\ 



( E - 2£i + fii y/UJh 



M 



£ - 2e N + n N ) 



(19) 



It is worth noticing that Richardson's basis, being of lower dimension, 
only yields the so-called seniority v = solutions. In contrast, our basis 
(|14"D provides the whole spectrum (any seniority) of the pairing hamiltonian 
eigenvalues: in the present case, of course, the v = 2 solutions are trivial, 
but when several pairs are present the v ^ solutions are important. 



In general the roots of Eq.(|l3|) cannot be given analytically for M > 4. 
Only the lowest collective eigenvalue, when G is large, turns out to be (in 
dimensional units) 



j-, G^oo 

£j r . ~ 



u=l 



-GQ 



(20) 



for any potential well. 

The other Af — 1 solutions, as is well-known, get instead "trapped" be- 
tween the unperturbed energies 



2e v -i<E v <2e v v = 2, •••TV (21) 

tAn alternative option uses for the basis the unnormalized vectors <£>'„ = In this 

case the A/"-dimensional matrix representing TL is filled with ones everywhere, except in 
the principal diagonal. 
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and in the limit G — > oo are given by the zeros of f(E/G) and depend upon 
the Af single particle energies e u and their degeneracies fl u . 

Remarkably, when one degeneracy, say Q u , becomes very large, then the 
"trapped" eigenvalues, for any value of G, coincide with the free ones, i.e. 



lim E„ 



2e^_! for 2 < ji < v 
2e, 



(22) 



n v ^oo >~ l 2e^ for v < \i < jV , 

whereas the lowest collective energy tends to — oo. 

Concerning the eigenf unctions, those (normalized) corresponding to the 
eigenvalues £ v = 2e v {y — 1, • • • AT) (seniority v = 2 states) in the basis (|il]) 
are 



1 



1 



$(°) 



(23) 



and describe a free pair sitting in the level j u . 

The f = eigenfunctions are more conveniently expressed in the basis 
Here they read 



Af 



^=o($*) = E y/^uWuK 



the coefficients tw„ fulfilling the system of equations 



Af 



(£ - 2e u ) w v + E = • 

The above is easily solved and yields the noticeable formula 

£-2e x 



£-2e l 



Since (P6|) entails for large £ 



(24) 



(25) 



(26) 



(27) 



the collective eigenstate in the G — » oo limit corresponds to a coherent 
superposition of all the s-quasibosons reading 



U^c) G ~°° E = E E 

y=l m, y =l/2 



(28) 



u=l 
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which is completely symmetric in the exchange of any pair of tp, thus exhibit- 
ing the same symmetry of the pairing Hamiltonian. 

Also from fl2"6|) one sees that, in the limit G — > where S u ~ 2e u , only one 



component of the basis, i.e. the z/-th one, survives in the wavefunction of the 
"trapped" states. The same occurs when one degeneracy, say Q u , becomes 
very large: as for the eigenvalues, the eigenvectors then coincide with the 
/x-th component of the basis for \i > v and with the (fi — l)-th component 
for fi < v. 

We now choose the harmonic oscillator potential for the single particle 
energies 

e N = (n + D %uj N = 0, • • • , oo (29) 
and for the degeneracies of the states available to a pair 

n N = ^(N + l)(N + 2) . (30) 

In this instance the secular equation (|T3"D, when the lowest Af levels are 
considered, becomes 

^ (JV+l)(iV + 2) _ 1 (31) 
^ Q 2N + 3-E G 

where G = G/2hu and the energies are measured in units of Tiuo {2e^ = 
2N + 3). 

In Fig.l the numerical solutions of ([31]) are displayed for J\f=6 and 8 ver- 



sus G. Remarkably, the dependence upon G is lost for G > 0.1. Furthermore 
in this regime the number of trapped solutions is obviously fixed by Af, but 
their dependence upon Af is quite mild. 

We now conjecture the Af trapped solutions of fl3~lD , which we label with 
the index N, to depend parabolically upon the latter, namely 

E N = aN 2 + bN + c N = 0, • • • AT - 2 (32) 

(the collective solution E c will be separately treated). 



To fix the coefficients a, b and c, we recast ( |31|) in the polynomial form 
E M + a{E M ~ x + a 2 E M ~ 2 + ■■■ a^E + a M = (33) 
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Figure 1: The figure shows the solutions E of Eq. fl3l|) , for 7V=6 (upper 
curves) and 8 (lower curves), as functions of G. One can see that with 
the harmonic oscillator well each trapped solution for G > 0.1 tends 
approximatively to the single particle energy 2e^. 

and compute the first three coefficients. They turn out to be 

a x = ~Af(Af+2)[G(Af+l)-3] (34) 
o 

a 2 = ^Af{Af- l)[—G(Af+ l){Af+ 2)(2A/" + 3) + 3J\f 2 + 1UV + 11] (35) 

a 3 = ^Ar(Ar-l)(Ar 2 -4)[G(Ar+l)(15Ar 2 + 40AT + 27)-15(AT 2 + 3Ar + 3)] . 

(36) 

Then the first three Viete equations, namely 

-ax - E c (37) 
a\ - 2a 2 - £ 2 (38) 
— 3a 3 — ai(a 2 — 3a 2 ) — E? c , (39) 



Af-2 

7V=0 
Af-2 

E^l 

7V=0 
AT- 2 

N=0 



s 



yield a non-linear system in the unknowns a, b and c, if E c is known. This 
system can be solved by expressing, via eq . (|37D , c as a function of a and b 



c(a,b) = --^h c + b -(Af-l)(Af-2) + ^(Af-l)(Af-2)(2Af-3) 
+ ^M(Af+2) [G(AT+1)-3]| . (40) 
In turn (50), inserted into (BB), yields b as a function of a. One finds 



15a(A/" 4 - 6A/" 3 + 13A/" 2 - 127V" + 4) ± VA 
6(a) = 15CAT- 1)^-2) (41) 



with 



A = _i5(Ar-l) 2 (AT-2){a 2 (Ar-l) 2 (Ar-l)(Ar-2)(Ar-3) (42) 



+ 20 



9E 2 C + QE C {N + 2){GM + G - 3) - 3(A/" 3 - 4A/" 2 - 13A/" - 11) 



- G 2 M(M - 2) (M + l) 2 (M + 2f + 3G(AT + 1) (AT + 2) (AT 2 - AM - 3)] } 

Finally, from (|39|) , an equation for a follows, not reported here, being 
quite cumbersome. While we have analytically solved (|37|) and (|38|) (they 
are of first and second degree in c and b, respectively), the non linear equation 
( |39| ) for a can only be solved numerically. Note that in ( |4TD the plus sign 
in front of the square root should be taken: however both options, when 
inserted in ([39]), lead to the same equation for a. 

It should be stressed that, because of the high degree of non-linearity of 
the latter, its solutions, when G is large, turn out to be extremely sensitive 
to the collective energy E c . To illustrate this issue we display in Fig. 2, for 
Af = 5 and G = 5, the surface (see (|39l)) 

S(a, E c ) =^E 3 R + 3a 3 + ax{a\ - 3a 2 ) + E\ (43) 

N=0 



whose zeros clearly give the values of the parameter a entering into (|32[) . It 
appears from the figure that even a tiny variation of E c induces a gigantic 
variation in S, implying that E c should be fixed with extreme precision in 
order to obtain the correct results for a, b and c. 
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-340 



Figure 2: The figure shows the surface S(a,E c ) for —342 < E c < —340 and 
—5 < a < 5. Here the number of levels is M = 5 and the pairing coupling 
constant is G = 5. 



Now a very good expression for the collective energy, obtained by expand- 
ing Eq. ([H]) in the parameter (2N + 3) / E and retaining the leading order, 
reads 

E^ = ~N{N +l)(N + 2) + ~(N + 1) . (44) 



Table 1 demonstrates the validity of ([44]) . Yet, its level of precision is 
not sufficient to allow a reliable determination, through (|39|) , of a, because of 
the dramatic non-linearity displayed in Fig.2. We are thus forced to improve 



upon (|44]) . To this purpose we insert into ([31]) the following expression for 
the collective energy 

E c = E^ + § (45) 
and again expand in the very small parameter 5/M(N) where 

M(N) = 2N + 3- Ej? ] . (46) 
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G = 1 


G = 5 


N 




£(0) 




E (e) 


E (0) 


E (D 


2 


-3.6055513 


-3.5 


-3.6041143 


-35.5192213 


-35.5 


-35.5192121 


3 


-14.096330 


-14.0 


-14.095861 


-94.0182425 


-94. 


-94.0182392 


4 


-32.582181 


-32.5 


-32.582012 


-192.515882 


-192.5 


-192.515881 


5 


-61.070599 


-61.0 


-61.070528 


-341.013793 


-341. 


-341.013793 


6 


-101.56156 


-101.5 


-101.56152 


-549.512104 


-549.5 


-549.512104 


7 


-156.05444 


-156.0 


-156.05442 


-828.010748 


-828. 


-828.010748 


8 


-226.54874 


-226.5 


-226.54873 


-1186.50965 


-1186.5 


-1186.50965 



Table 1: Comparison between the exact (e) and the approximate (0) (eq. (f44"D ) 
and (1) (eq.([47])) collective energies for some values of N and G=l and 5. 



We thus obtain, in lieu of 



the expression 



E m 



E^ + 



JV-l 



± v (N + l)(N + 2) 
G 2iV + 3-^ 0) 



^ x (jV + l)(AT + 2) 
n=o (2N + 3 - E { c ] ) 2 



-i 



C{M, Ef>) + 2Af (E® + M + 2) - 8/G 



(OK 



(47) 



8E" 



where 

C(Af,E^) 



1 



3-Ei°y 



'3 -Ei°y 



(48) 



(\l/ being the Digamma function). Eq.([|7]) provides the exact collective en- 
ergy, for all practical purposes, as shown in the fourth and last column of 
Table 1. 

However, for weaker G, (|47|) becomes less accurate, because the expansion 
parameter (2iV + 3)/E is no longer small. But this occurrence is of no 
consequence for the trapped solutions (up to about G = 1), since a weaker G 
also means a less severe non linearity. If, however, a more accurate value for 
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the collective energy is wished, it can be found through the recursion relation 



1 

G 



^ (N + l)(N + 2) 
2N + 3-E { c k) 



N=0 



^ (N + l)(N + 2) 
n=o (2N + 3-E { c h) y 



(49) 

the index k labelling the order of the iteration. We have numerically checked 
that the iterative expansion ( f49|) converges to the exact solution for G > 0.2. 

However, the physical interesting domain for the pairing problem occurs 
for (5 = 0.05-0.1, since [| 



G 



27 1 A MeV for protons 
22 /A MeV for neutrons 



(50) 



A being the nuclear mass number, and hu should correspond to the aver- 
age distance between the single particle levels inside the last occupied shell, 
namely « 0.69 MeV in Pb and « 0.74 MeV in Sn. 

Interestingly an accurate expression for the collective energy can also be 
given in the regime < G < 0.2. Indeed here of E c (G,J\f) we know the value 
for G = (namely 3) and where it vanishes, namely for 



Go 



^(N+1)(N 

.N=0 



-, -1 



(2N + 3) 

Moreover, the value of its derivative is —2 in G = and 
dE 



(51) 



dG 



'■^ 1 (iV + l)(V + 2) 

.N=0 



2N + 3 



2 r 



^ (N + l)(N + 2) 
L ko (2iV + 3) 2 



-, -i 



(52) 



in G = Go. 

These four constraints are fulfilled by the cubic 



£?(°)(G) = 3-2G-^9 



dE, 



dG 



G 



6 + 



dG 



G 



(53) 

which thus provides an excellent starting point for the evaluation of the 
collective energy. Proceeding indeed as done in the large G domain, a per- 
turbative expansion in a parameter 5 can be set up, leading again to formula 
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G = 


0.05 






M 


E (0) 




E W 


E f) 


E W 


E W 


2 


2.87617 


2.88394 


2.88349 


2.88348 


2.88348 


2.88348 


3 


2.81180 


2.87628 


2.86196 


2.86014 


2.86012 


2.86012 


4 


2.70198 


2.89185 


2.84845 


2.82471 


2.82056 


2.82046 


5 


2.53580 


2.86975 


2.80333 


2.75484 


2.74104 


2.74028 


6 


2.24594 


2.63571 


2.54951 


2.52958 


2.52886 


2.52886 


7 


1.61587 


1.84837 


1.83108 


1.83094 


1.83094 


1.83094 


8 


.134714 


.136135 


.136135 


.136135 


.136135 


.136135 








G = 


0.1 






M 




E W 


E f) 


E f) 


E W 


EW 


2 


2.70757 


2.72972 


2.72822 


2.72822 


2.72822 


2.72822 


3 


2.45586 


2.61626 


2.58336 


2.58335 


2.58335 


2.58335 


4 


1.96617 


2.21918 


2.18238 


2.18238 


2.18238 


2.18238 


5 


0.919942 


1.00125 


1.00000 


1.00000 


1.00000 


1.00000 


6 


-1.66966 


-1.47625 


-1.48025 


-1.48025 


-1.48025 


-1.48025 


7 


-8.37559 


-4.88888 


-5.44566 


-5.44568 


-5.44568 


-5.44568 


8 


-24.4931 


-3.26859 


-10.7011 


-11.0491 


-11.0546 


-11.0546 



Table 2: Comparison between the exact and the approximate 
[eq.(|5^) and eq.([49"|)1 collective energies for some values of Af and G=0.05 
and 0.1 



(P), but with the input E^ now given by the k-th iteration of (|53"D . How 
remarkably accurate the results we obtain are, is displayed in Table 2. 

With the collective energy fixed, the coefficients a, b and c can be found. 

It should be reminded that, being the system of the equations (0), (|38|) 
and (^) non-linear, more than one set of solutions is generally found: how- 
ever the appropriate set is easily selected, being the one yielding energies in 
between the single particle levels. 

We quote in Table 3 our predictions for the eigenvalues of the pairing 
hamiltonian for one pair in the M = 5 case, using as input (0) when G = 1 
and G = 5 and ( |53D when G = 0.05 and G = 0.1. These leading orders 
are iterated via (E^) until self-consistency is reached. Our results are seen 
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G = 


0.05 


G = 


0.1 


G -- 


= 1 


G -- 


= 5 


N 


E ( - e) 

N 


7?(<W) 


N 




N 




E { - e) 

N 







4.2872 


4.2812 


3.4245 


3.4266 


3.1583 


3.1601 


3.1493 


3.1510 


1 


6.0892 


6.1056 


5.6302 


5.6237 


5.3673 


5.3621 


5.3524 


5.3472 


2 


8.1171 


8.1021 


7.8422 


7.8485 


7.6136 


7.6190 


7.5965 


7.6015 


3 


10.266 


10.271 


10.103 


10.101 


9.9314 


9.9297 


9.9157 


9.9140 



Table 3: Comparison between exact (e) "trapped" solutions of Eq.(|31|) and 
approximate (app) ones, obtained from the ansatz (^) for M = 5 levels. 
The coefficients (a, b, c) of the parabola are (0.086,1.738,4.281) when G = 
0.05, (0.014,2.183,3.427) when G = 0.1, (0.027,2.175,3.160) when (5 = 1 and 
(0.029,2.167,3.151) when G = 5. 



to agree with the exact ones obtained via the numerical solution of ( pTj ) to 
better than 0.27%. 

In this letter, to pave the way to the problem of any number of fermions 
pairs, we have solved, almost analytically, the pairing problem for one pair 
living in any number of harmonic oscillator levels. Crucial for this achieve- 
ment has been the conjecture Q32"D , possibly related to the specific degeneracy 
of the harmonic oscillator single particle levels. If this is true formulas like 
( |32"D can as well hold valid for others one body potentials, providing their 
degeneracy is known. 
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